Landau-Zener dynamics of a nanoresonator containing a tunneling spin 
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We study the Landau-Zener dynamics of a tunneling spin coupled to a torsional resonator. For 
strong spin-phonon coupling, when the oscillator frequency is large compared to the tunnel splitting, 
the system exhibits multiple Landau-Zener transitions. Entanglement of spin and mechanical an- 
gular momentum results in abrupt changes of oscillator dynamics which coincide in time with spin 
transitions. We show that a large number of spins on a single oscillator coupled only through the 
in-phase phonon field behaves as a single large spin, greatly enhancing the spin-phonon coupling. 
We compare purely quantum and semiclassical dynamics of the system and discuss their experimen- 
tal realizations. An experiment is proposed in which the field sweep is used to read out the exact 
quantum state of the mechanical resonator. 
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I. INTRODUCTION 

The Landau-Zener modeP describes a two-state sys- 
tem in which the bias between diagonal states varies lin- 
early with time as they are swept through an avoided 
crossing. It is one of the few practically important time- 
dependent Hamiltonians for which the Scrhodinger equa- 
tion is exactly solvable. The Landau-Zener method has 
recently found a natural application in the experimental 
characterization of single molecule magnets^ . Theoretical 
studies of Landau-Zener transitions in nanomagnets have 
included many-body effect^ and superradianc^. Some 
important theorems have been proven about generaliza- 
tions of the Landau-Zener probleirP, and certain multi- 
level cases have been exactly solvecP. It has been used as 
a model for the dynamics of quantum phase transition^, 
and topological defect formatioiP. 

A natural extension of the two-level quantum physics 
is the two-level system coupled to one or several quan- 
tized modes of a harmonic oscillator. Studies of Landau- 
Zener os cillato r dynamics have probed coherent^ 
dissipativ^Sm^ and temperature-dependentP^ effects. 
Landau-Zener interferometry has provided a quantitative 
measure of coupled dynamics— and has been experimen- 
tally verified in the nanomechanical measurement of a 
superconducting qubitP^. 

The Landau-Zener effect in spin systems should be 
considered in conjunction with the transfer of angular 
momentum manifested in the Einstein - de Haas effect. 
This effect has allowed precision measurement of the 
magneto-mechanical ratio of a thin ferromagnetic film on 
a microcantileveili^. Torsional oscillators have been used 
as precision torque magnetometers in nanomechanical de- 
tection of itinerant electron spin-flip at a ferromagnet- 
normal metal junctiorP^and measurement of phase tran- 
sitions of small magnetic disks in and out of the vortex 
statd^. Semiclassical models of Landau-Zener dynam- 
ics have been developed to describe magnetic molecules 
coupled to mechanical resonators and bridged between 
conducting lead^^^Hni, A full quantum treatment of the 



interaction between a single sp in and a torsional oscilla- 
tor has recently been developecP^^. 

Realizing a quantum magneto-mechanical system with 
strong spin-phonon coupling has been an experimental 
challenge. A recent experimenll^^] has shown the first evi- 
dence of strong spin-phonon coupling in a single molecule 
magnet grafted onto a carbon nanotube. Spin reversal of 
the single molecule magnet during a Landau-Zener sweep 
coincides with an abrupt increase in the differential con- 
ductance through the carbon nanotube. This has been 
interpreted as the spin transition exciting a longitudi- 
nal stretching mode of the carbon nanotube, which en- 
hances electron tunneling from the lead onto the nan- 
otube through electron-phonon coupling. 

We propose multiple schemes to realize strongly cou- 
pled dynamics of a tunneling macrospin with torsional 
oscillations of a nanoresonator in a Landau-Zener ex- 
periment. We investigate the Landau-Zener dynamics 
of a tunneling spin coupled to a torsional oscillator, us- 
ing a fully quantum mechanical model. The oscillator 
could be a torsional paddle resonator, a microcantilever, 
a carbon nanotube, or a single magnetic molecule be- 
tween two point contacts. The tunneling spin could be a 
single molecule magnet, an ensemble of single molecule 
magnets, or a single-domain ferromagnetic particle with 
strong uniaxial anisotropy. For a collection of single 
molecule magnets placed on a torsional resonator or can- 
tilever far apart from each other that they are not directly 
coupled through dipole interactions, we develop a semi- 
classical model of magnetization dynamics. We predict 
superradiant enhancemenlPl of the spin-phonon coupling 
for this ensemble system. Comparison of these two mod- 
els shows their correspondence. 

The coupling between spin and mechanical angular 
momentum is mandated by the conservation of total an- 
gular momentum J = S + L, with L bcingthe mechani- 
cal angular momentum. In a free particle, when a spin 
tunnels from S to — S, the particle must change its me- 
chanical angular momentum L. This changes its kinetic 
energy by an amount of order h?S^/I, where Iz is the 
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moment of inertia about the rotation axis. For a macro- 
scopically large body, the large moment of inertia makes 
this rotational kinetic energy negligibly small. But for a 
small particle this can become comparable to the energy 
gain A due to tunnel splitting^^. The ratio of these two 
quantities, the magneto-mechanical ratio a = 21%^ / 
determines the ground state of the system for a free 
particle For large particles a <C 1 and the ground state 
is the well-known tunnel split state ^ ~ \ips) + IV'-s)- 
For small particles, such as the Feg single molecule mag- 
net with Iz ~ 10~^^ kg-m^, a ^ 1 and spin tunneling is 
suppressed as the spin localizes in either direction along 
the easy axis. 

Similar effects arise in systems that undergo torsional 
oscillations. Examples are a single molecule magnet 
bridged between conducting leads, a nanomagnet at- 
tached to a carbon nanotube bridge, or a nanomagnet 
coupled to a resonator such as a torsional paddle os- 
cillator or microcantilever. The mechanical resonance 
occurs at a frequency ujr = y/kJTz, where k is the ef- 
fective stiffness against the linear restoring torque and 
Iz is the moment of inertia of the nanomagnet-resonator 
combination. A convenient measure of the effect of os- 
cillations is the dimensionless parameter r = Scj^/A, the 
ratio between the oscillator energy to tunnel splitting. 
As we will see in Sec. the coupling between magne- 
tization and oscillator dynamics is given by the factor 
A ~ \/oJr — yj2hS'^ / IzUJr- The most interesting effects 
occur for strong coupling A ~ 1 and oscillator frequency 
much larger than tunnel splitting r ^ 1. A large spin, 
small moment of inertia, and weak torsional spring con- 
stant are required for strong coupling, article is organized 



II. MODEL 

A. Landau-Zener Transitions in a Two-State 
System 

We review relevant features of the Landau-Zener 
model, which describes a two-level system driven by a 
classical field that varies linearly in time. The LZ Hamil- 
tonian is 



. vt l\ 

tlLZ = -^^z - yCTa: 



(1) 



in terms of Pauli matrices Uz and Ux = cr+ -\- <J-, where v 
is the sweep rate and A is the tunnel splitting. Diabatic 
states I t) and | \) are eigenstates of cr^ with diabatic 
energies E^:^{t) = zLvt/2, which are the linear functions 
in Fig. [2a[ We take the sweep rate v positive, so the 
positive (negative) sign corresponds to spin down (up). 
For nonzero A, the diabatic states are not eigenstates 
of the Hamiltonian. Diagonalizing H^z gives adiabatic 
energies 



(2) 



which are the upper and lower curves in Fig. [2a| with 
splitting A at t = 0. The corresponding adiabatic eigen- 
states 1-1-) and I—) are 
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|±) = ^(C^|t)TC±U)), 



(3) 



where C± depend explicitly on time. 




FIG. 1. Possible experimental geometries described by the 
models studied in this paper. In both cases the easy axis of the 
macrospin coincides with the rotation axis of the oscillator, 
(a) Single molecule magnet grafted on a carbon nanotube. (b) 
Ensemble of single molecule magnets on a nanocantilever. 

as follows. In Sec. |ll]we briefly review the Landau-Zener 
model, and construct the quantum mechanical model of 
a spin coupled to a torsional resonator with an external 
magnetic field that varies linearly in time. Sec. 
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tains numerical and analytical results of the fully quan- 
tum spin dynamics for a variety of parameter ranges. 
Oscillator dynamics are presented in Sec. |IV| A semi- 
classical model of superradiant dynamics in an ensemble 
of spins on a single resonator is developed in Sec. |V] 
Finally, we discuss the interpretation of our results for 
various experimental realizations in Sec. VI 
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ypF+A2' 
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For times \t\ ^ A/v the adiabatic states asymptotically 
coincide with the diabatic states. 
The state of the system 



(5) 



evolves according to the time-dependent Schrodinger 
equation 



ih 



~dt 



(6) 



with initial conditions c^{—oo) ~ 0, |c^(— oo)| = 1. After 
eliminating c^, we obtain the second order differential 
equation 
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ctit) = (7) 



which can be put into the standard form of the Weber 
equation. The exact solutiorP gives 



(i) = V^e-^^/42?_„_i(-iz) 



(8) 
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where 
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V = 27, 



(9) 



and -D_i/_i(— jz) are parabohc cyhnder functions. The 
staying probabihty for the spin-down state as function of 
time is P{t) ~ |c^(t)p. The exact asymptotic hmit for 
t = oo, known as the Landau-Zener probabihty, is 



LZ 



7^ 

2hv 



(10) 



P{t) and Plz are shown in Fig. 2b The same P{t) 
and Plz can be obtained from the Heisenberg equa- 
tions of motion for {az{t)). An intuitive understanding 
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FIG. 2. (a) Adiabatic E±{t) and diabatic Efi{t) energy levels 
of the LZ Hamiltonian as a function of time, (b) Probability 
P{t) of staying in the initial | \) state as a function of time, 
and asymptotic staying probability Plz- 

of the Landau-Zener transition comes from considering 
the time spent in the tunnehng region between adiabatic 
states and the tunneling time between these states. Let 
tlz ~ max{y/h/v , A/v) be the time spent in the tunnel- 
ing region and ta ~ H/ A be the tunneling time at the 
crossing. The Landau-Zener exponent is proportional to 
the ratio of these times e ~ tlz/ta- For a slow sweep the 
system will evolve adiabatically, spending long enough in 
the tunneling region that it will continually relax to the 
ground state, making e 1 and Plz — >■ 0. In the op- 
posite limit, a fast sweep through the tunneling region 
makes e ^ 1 and the staying probability saturates at 
Plz ^ 1. 



B. Landau-Zener Transitions in a Spin-Oscillator 
System 

Consider a tunneling spin which is projected onto the 
lowest tunneling doublet. This spin is coupled to a 
torsional nanoresonator with rigidity k that can rotate 
about the z-axis, which coincides with the easy axis of 
the spin. The Hamiltonian i j^^ l ^^l , 



- h^Ll Izl4£ W{t) A 

(11) 

The fundamental frequency of torsional oscillations is 
= y^k/Iz, where Iz is the moment of inertia of the 



resonator about its rotation axis. An external longitu- 
dinal magnetic field Bz{t) applied along this axis cre- 
ates a time-dependent energy bias W{t) — 2SgiXBBz{t)- 
The Landau-Zener problem describes a linear field sweep, 
W{t) = vt. The operator of mechanical angular mo- 
mentum, Lz = ~id^, and the angular displacement (p 
of the oscillator obey the usual commutation relation 
[(t),Lz\ = i. 

The last term in the Hamiltonian describes the entan- 
glement between spin transitions and mechanical rota- 
tions. A typical single molecule magnet has a large spin 
and strong uniaxial anisotropy, producing a zero-field 
splitting between degenerate ground states \'ip±s) point- 
ing in either direction along the easy axis. Any symme- 
try breaking interactions, such as transverse anisotropy 
or an external field, break this degeneracy producing tun- 
nel split states ^' ~ jV's) ± IV'-s) which are represented 
by the pseudospin cr. The tunnel splitting A is generally 
many orders of magnitude less than the energy to the 
next spin level. In the case of the spin- 10 single molecule 
magnet Feg, the crystal field Hamiltonian describing the 
magnetic anisotropy is Hs = —DS^ 
Full perturbation theorjESl gives 



dS^, with d< D. 



A ^ 



8S^ 
7ri/2 



_d_ 
AD 



(12) 



where we can see that A ^ 2SD, which is the distance 
to the next spin level. The crystal field Hamiltonian Hs 
is defined with respect to coordinate axes that are rigidly 
coupled to the molecule or crystal. Because the particle 
is free to rotate, the crystal field Hamiltonian must be 
transformed to the fixed frame of the laboratory. Pro- 
jecting the crystal field Hamiltonian onto the lowest tun- 
neling doublet, rotating to the lab frame using U{Sz) = 
e»s.0^ where Sz\^±s) ^ ±5|i^±s), H'g = UHslf-^ gives 
the final term of the Hamiltonian. 

We now consider the spin-oscillator Hamiltonian with 
a linear field sweep W{t) = vt. Introducing the usual 
annihilation and creation operators, a and a\ 

^=/^(«t+„), Lz=^^(a^-a) (13) 



into Eq. ( 11 ) gives 



(14) 

where we have dropped unessential constant terms. We 
will find it useful to adopt dimensionlcss units tl' — H / A 
and t' = At/h, 

H' = ra^a - "^az - ^(e-'^('^'+''V+ + e^^^-^'+^V.), 

(15) 

which shows that the system depends on three parame- 
ters. The parameters 



l2hS^ 

IzUJr ' 



huJr 



(16) 
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describe the spin-oscillator relationship. A is the coupling 
strength between the spin and oscillator and r is the ratio 
of mechanical oscillation to tunnel splitting frequency. 
The relationship between A and r can be understood by 
the so-called magneto-mechanical ratio, 



= A^r = 



2^ 



(17) 



which is the ratio of the change in rotational kinetic en- 
ergy associated with a spin transition S — S to the 
tunnel splitting energy. The third parameter is the ef- 
fective sweep rate w', or equivalently the Landau-Zener 
exponent e defined in Eq. (10), 



Ye 



hv 
A2- 



(18) 



We choose the spin up/down basis for the two-level 
system and a Fock state basis for the harmonic oscilla- 
tor. A direct product of these two bases will form the 
basis of the spin-oscillator system. Matrix elements of 
the Hamiltonian Eq. ( 14 1 are 



(19) 



A! 



where a — —1,1 corresponds to spin down and up states, 
respectively. The full Fock space has an infinite num- 
ber of states, although we will use a truncated basis for 
numerical computations. Tunneling matrix elements 



A. 



A/ 



.(A), 



(20) 



depend on the coupling A through matrix elements of the 
displacement operator D{i) = exp(^a^ — i = —iX, 

^mn{X) = e-^'/2(-*A)™-'s/44"""'(A') (21) 

V m! 

for m > n, and m ^ n for m < n. L^™ (x) are gener- 
alized Laguerre polynomials, and the real parameter A is 
defined in Eq. (16 1. The first few are 



Kii 



(l_A2)e-^V2. 



(22) 



III. LANDAU-ZENER SPIN-OSCILLATOR 
DYNAMICS 

A. Adiabatic energy levels 

Numerically solving det{H — EI) = gives the adia- 
batic energy levels En±, shown in Fig.jSj Diabatic energy 



levels Eni^, dotted lines in the insets of Fig.lSj are eigen- 
values of the noninteracting part of the HamTTtonian (the 
first two terms in Eq. (14)), given by 



E., 



nr ± 



v't' 



(23) 



A 2 

The spin down (up) states have positive (negative) slopes 
with y-intercepts nuj. Diabatic energies E„i and Em-]- 
cross at times 

t' • 



k — , k^m-neZ. (24) 
v' 



When the oscillator frequency is much larger than the 
50r 

40 : 



til 




t' 



FIG. 3. Energy (in units of A) as a function of time for 
r = 20, A = 1, e = 1.35. Solid lines are adiabatic energy 
levels En±, and diabatic energies Enif are dashed lines in the 
insets. Crossings occur at tk- 

sweep rate, r ^ v' , the transitions are independent. Note 
that the indices on the adiabatic and diabatic energies 
only coincide near t = 0, but will in general be different 
after successive crossings. The tunnel splittings |Amn| 
between adiabatic states occur at the crossing of dia- 
batic energies E^i and Enu and depend on the coupling 
strength through Eqs. (201 and (21 1. When r > v' , 



successive transitions occur within short times of each 
other. Once r < v' there are many closely spaced levels 
near t = Q. 

Consider a single spin initially spin-down with the os- 
cillator in the zero phonon state, i.e. ^(< = — oo) = |0)| | 
) . The system is initially in the adiabatic energy state 
_Eo- which corresponds to the diabatic state Eq^. At 
to = 0, diabatic states Eq^ and Eq^ cross, and adiabatic 
states Eq- and Eq+ approach each other with minimum 
separation |Aoo| — Ae^"^ If the spin remains in the 
initial adiabatic state Eq- after the avoided crossing, it 
flips and will see no more possible transitions, as Eq_ 
coincides with ii^ot long after the avoided crossing at t^. 
If the spin does not flip, it will follow the adiabatic state 
Eq+ which coincides with Eq^ long after to- The next 
crossing between diabatic states Eq^ and Ei^ occurs at 
ii, with tunnel splitting |Aoi| = Ae"'*' ^^A between di- 
abatic states Eqj^ and . Remaining in the adiabatic 
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state Eq+ will coincide with Ei^ for times long after ti . If 
the spin does not flip at ti, the system will remain in the 
adiabatic state, coinciding with Eq^ long after ti. 
In general the crossing between state |0)| \.) and f) 
occurs at tk with splitting |Aofc| — Ae~^ ''^|'«ofe(A)|. No- 
tice that the avoided crossing between Ei_ and Ei^ at 
to — 0, given by |Aii| = Ae^'^ /^|1 — does exactly go 
to zero when A = 1. 



B. Strong coupling 

We study the dynamics of the spin-oscillator system for 
various parameter ranges. Expanding the wave function 
of the system in this basis 



(25) 



m=0 cr=±l 



the time-dependent Schrodinger equation yields the sys- 
tem of coupled differential equations, 



.dC„ia , v't' 



(26) 



E 



We solve this system of equations numerically with a 
truncated oscillator basis. First we consider the initial 
state of the spin system to be spin-down with the oscilla- 
tor in its quantum ground state |^'(— oo)) — |0)| |), which 
gives Co__i(— cxd) = 1 with all other Cm.a{—oo) = 0. 

Strong coupling (A ^ 1) of spin dynamics to torsional 
oscillations results in rich dynamics of both the spin and 
the oscillator. Calculating the expectation value of cr^. 



(27) 



we define the probability of staying in the initial spin- 
down state as 



(28) 



A comparison of staying probabilities for different param- 
eters is shown in Fig. [4] For r 3> 1, the spin transitions 
are clearly independent, as shown in Figs.[4a|and|4b[ The 
tunnel splitting at each crossing is strongly renormalized, 
according to Eq. (201, which leads to strong dependence 



of the transition probability on the coupling. 

Consider the crossing of diabatic energies Emi and 
En^. For the system initially in the |r7i)| |) state, which 
corresponds to the lower of the two adiabatic states long 
before the avoided crossing, the probability that the sys- 
tem will stay in the initial state is 



Prr 



7rA^|Km»P 
2hv 



(29) 
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FIG. 4. Time dependence of the probability that the spin 
stays in the initial spin-down state for initial state >!'(— oo) = 
|0)| \) with e = 1.35. Vertical lines at tk denote avoided 
crossing of adiabatic energy levels. Horizontal lines are exact 
results Pm for independent transitions. 



When the system is initially in the |0)| 1) state, all dia- 
batic crossings will occur between energies E^i and En-\- 
The transition probability Pon — e~'^°" at each crossing 
depends on |KonP- Using Lq{x) = 1 we obtain 



ttA^-^ A2" 



(30) 



2hv nl 

After the first avoided crossing at to = 0, the asymptotic 



6 



staying probability in the initial state is Pqo = e"*^"". The 
next avoided crossing occurs at ti , and the probability of 
staying in the spin down state after ti is Pqi = e^*^"!^ 
Thus the total staying probability after two avoided 
crossings is Poo^bi- We define P/v as the probability of 
remaining in the initial state after N avoided crossings, 

Pjv = exp - ^ eo« I • (31) 

\ n=0 I 

0, we recover the exact Landau Zener 



In the limit N H 
probabihty P^z, 



lim p 

Af->-oo 



7rA~e 



N 



exp 



2hv 




exp 



2hv 
(32) 

Fig. 4b shows staying probability for larger coupling, 
A = 2. We see that as the tunnel splitting at each 
avoided crossing is more strongly renormalized, it takes 
more crossings to reach the final Landau-Zener probabil- 
ity. 

As the oscillator frequency decreases compared to the 
sweep rate, r > 1, the transitions are no longer com- 
pletely independent, although small oscillations about 
individual plateaus can still be seen in P(i). This is be- 
cause the transitions happen within a small multiple of 
the Landau-Zener tunneling time tj^z- When the oscilla- 
tor frequency and tunnel splitting are close to resonance 
r 1, the transition probability initially approaches Plz 
and then shows collapse and revival behavior around this 
limit, as shown in Fig. [4cj For r <C 1 the revivals become 
much weaker and the probability resembles the tradi- 
tional LZ probability. 
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FIG. 5. Time dependence of the probability that the spin 
stays in the initial spin-down state for an initial coherent os- 
cillator state ^(— oo) = |/?)| \) with e = 1.35. Vertical lines 
at tk denote avoided crossing of adiabatic energy levels. Hor- 
izontal lines are exact results Pjv for independent transitions 
when starting in the |0)| \) state. 

When the oscillator is initially in a coherent state |/3), 

^{-^)^mi)=e-\''\ (33) 



where the complex number /3 = |/3|e'* is proportional to 
the amplitude of initial oscillations. When /3 ^ 1 the 
spin transitions follow approximately the same asymp- 
totic values Pat as the quantum ground state case. For 
/3 < 1 the staying probabilities, an example of which is 
shown in Fig. [5] depend on the magnitude and phase of 
the initial coherent state. The maximum angular dis- 
placement and velocity of a coherent state are related to 
/3 through 



'/'max = 2A|/3|, 



dip 
df 



2rA|/?|, (34) 



where ip — 2S{<f)). 



C. Weak coupling 

When the spin dynamics of the nanomagnet are weakly 
coupled to its rotational dynamics A <C 1, there is lit- 
tle observable effect of rotations on spin flip probability. 
The first crossing that occurs at to = has tunnel split- 
ting Ago — Ae"'*' /^, which tends to unity for small A. 
When r ^ 1 the first transition at io = approaches 
^00 — e^'^'"'. The second independent transition occurs 
at ti approaches Plz^ although the difference between 
Poo and Plz is very small. When r ~ 1 the adiabatic 
transitions are no longer independent, and occur close to 
the Landau-Zener tunneling time interval. 



IV. OSCILLATOR DYNAMICS 

We compute the expectation value of the torsional ro- 
tation angle as a function of time 



A (C'm-i-l -xCn 



Vm + 1 + C*-^^i ,jCrn.aVrn) . 

(35) 

For strong coupling A ~ 1, the dynamics of the resonator 
shows a delay before the onset of large oscillations for 
r 1, which occurs at ti, shown in Fig. |6a[ When the 
coupling is stronger, Fig. |6b| shows many more changes 
in the oscillatory motion, consistent with more avoided 
crossings. The oscillation amplitude changes slightly at 
subsequent tk- As r decreases towards 1, the interval of 
large oscillations becomes shorter. When r < 1, there is 
a single transition region which gives way to harmonic 



oscillations, shown in Fig. 6c Near r — 1, the amplitude 
of oscillations tends to increase as r decreases for fixed 
A. 

The angular displacement of the torsional resonator 
also shows interesting effects even for small coupling. 
When r ^ 1, large torsional oscillations do not begin 
at the first crossing. This can be understood as follows. 
The t = crossing occurs between spin up and down 
states, both of which correspond to the ground state of 
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FIG. 6. Time dependence of the rotation angle expectation 
value for initial state ']/(— oo) = |0)| ^) with e = 1.35. Vertical 
lines at tk denote avoided crossing of adiabatic energy levels. 



the resonator. Although there is a small increase in dis- 
placement angle at this crossing, the largest increase oc- 
curs at the second crossing between |0)| i) and |1)| t) 
at ti. This delay agrees exactly with the semiclassical 
treatment by Jaafar et al.l^. Following this time the os- 
cillator is in a superposition of ground and excited states. 
When r <1 successive transitions occur in a short dura- 
tion compared to r^z, and there is no observable delay 
in the onset of oscillations. The oscillation amplitude de- 
pends on the sweep rate. Numerical results suggest that 
the largest amplitude oscillation reaches a maximum near 
e ~ 2 for A = 1 and r = 20. 

When the oscillator is initially in a coherent state 
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FIG. 7. Time dependence of the rotation angle expectation 
value for initial coherent state VE'(— oo) = |/3)| \) with e = 1.35. 
Vertical lines at denote avoided crossing of adiabatic energy 
levels. 



normal oscillations with maximum amplitude <y9niax = 
2A|/3| occur up to t-i, as shown in Fig. [7j At t„i the 
amplitude decreases slightly and decreases again at to- 
A large increase occurs at ti, similar to the case where 
the oscillator is initially in its quantum ground state. The 
amplitude of oscillations after ti tends to be larger when 
the oscillator is initially in a coherent state, but not by a 
large amount. We observe a subsequent change in oscilla- 
tion amplitude at t2 and ^3. The oscillator dynamics are 
not as sensitive to the initial phase of the coherent state 
as the spin dynamics, although there is some variation in 
maximum amplitude. 



V. COLLECTIVE DYNAMICS OF SPINS 
COUPLED TO A MECHANICAL RESONATOR 

Consider, instead of a single nanomagnet, an array of 
single molecule magnets with their easy axes mutually 
aligned with the axis of rotation of the resonator. If they 
are far enough apart that dipolar coupling is negligible, 
they will only be coupled through the effective field due to 
torsional oscillations. Because the angular displacement 
is the same for each molecule, this results in collective 
coherent dynamics, described by a variant of the Dicke 
Hamiltonian. For N single molecule magnets, we define 
the operator of total low-energy dynamics as 

A ^ 

Hr^--R,, R = ^t' (36) 

i=l 

where the index i labels each magnetic particle. Again 
transforming to the lab frame by performing a rotation 
by angle (f) to the lab frame, but now using the total spin. 
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we obtain 

= - ^ (cos (2S(j)) + sin (25?!)) i?^;) . (37) 

The full Hamiltonian for the array of single molecule 
magnets is 

^ ^h^Ll h^lcp^ W{t) 

2h 2 2 ^ 

- ^ (cos (250) + sin {2S4>) Ry) , (38) 
where W{t) — vt. The Hamiltonian can be written as 



HsR = Hose — 2-^<=ff ' 



(39) 



where i/osc is the uncoupled oscillator Hamiltonian and 



H 



eff 



5H 

Jr. 



A cos {2S(l))ex + A sin {2S(f>)ey + We^ 



(40) 

is the effective magnetic field. Noticing that Hsr is lin- 
ear in Rx,Ry,Rz, we can see that [II'^,Hsr] = 0, so 
= R{R + 1) is a conserved quantum number and R 
behaves as a single large isospin. We are interested in 
the maximum value of R, i?max — N/2, which can be 
experimentally realized by preparing the system with a 
strong longitudinal magnetic field such that all spins are 
pointing down. 

The Heisenberg equations of motion ihdA/dt — [A, H] 
are 



'hLu^4) - A5 (sin {2S(i))R^ - cos (2S(j))Ry 



HR^ = WRy - A sin {2S(I))R,^ 

hky = -WR^ + Acos i2S(j))R, 

hRz = -A cos {2S(f))Ry + A sin {2S(I))R^. 



(41) 
(42) 

(43) 
(44) 
(45) 



discussed in Sec. H A Second, these equations are not 
independent, but 



— R2 = 
dt 



(47) 



which is equivalent to R^ = constant, which we had 
found as a constant of motion of the Hamiltonian. Be- 
cause the length of R is fixed and large in magnitude, 
we see that the equations of motion for R are equiva- 
lent to the Landau-Lifshitz equations for a classical spin 
of fixed length precessing in a magnetic field. Dividing 
Eqs.(41)-(45) by R shows that the direction of the to- 



tal spin follows the Landau-Lifshitz equation, which is 
mathematically equivalent to the Schrodinger equation 
of a spin-half particle precessing in a magnetic field, 



h— cr X Hcff , 
dt 



R 



(48) 



The equations of motion for Rx,y,z can be divided 
through by R to giving identical equations of motion for 
a pseudospin cr = R/R of unit length. Substituting this 
into the equation of motion for (j) and eliminating gives 
a second order equation of motion for for the dynamics 
of the resonator. 



+ r'^ip = ~aR (sin {tp)ax — cos i(p)cry) , (49) 



where ip — 2S(j), the prime denotes derivative with re- 
spect to dimensionless time t' = At/h, r and a are de- 
fined in Eqs. (16) and (17 1, respectively. 

( [49| shows that the spins ex- 
This is a simple 



The right hand side of Eq. 
ert a collective torque on the resonator, 
yet meaningful result. The equation of motion is simi- 
lar to the semiclassical treatment of a single spirPSl, but 
with the torque on the resonator increased by a factor of 
R. Because the amplitude of oscillation is proportional 
to the number of magnetic molecules N , this can be in- 
terpreted as a signature of Dicke phonon superradiancd^. 
For a simple harmonic torsional oscillator, the phonon 
field is the angle of displacement from equilibrium 0, and 
the driving torque is proportional to i? = N /2. 

Returning to the quantum model we see that for the 
case of superradiance, a aR, and A — \J a/r becomes 



The equations of motion show a few important proper- 
ties. First, the time derivative of the z-component of the 
total angular momentum equals the elastic torque 



d 
dt 



(46) 



If the spin-rotor system were completely uncoupled from 
its environment the total angular momentum, spin plus 
rotational, would be conserved. In the limit 0, we 
would obtain Heisenberg equations of motion for Rx,y,z- 
Solving this system of equations gives the same Landau- 
Zener probability of spin flip as the Schrodinger picture. 



A 



SR 



'R\(x VNX 



(50) 



This provides a viable method of increasing the coupling 
in a realistic experiment, by increasing the number of 
individual nanomagnets on the resonator. The usual dif- 
ficulty of realizing strong coupling is that reducing the 
moment of inertia by even two orders of magnitude has 
a small effect on the coupling due to the inverse quar- 
tic root dependence of the coupling on the moment of 
inertia. 



The Heisenberg equations of motion, Eqs. (41 )-(45 1 are 



operator equations which should be averaged over the 
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quantum state of the system. Since the spin R is classical 
the averages decouple, such as {si'D.(ip)ax) {siTi{ip)) (ax) 
in Eq. (49), which yields classical-like equations of mo- 
tion. We solve these equations of motion numerically. 
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FIG. 8. Time dependence of the effective probability P{t) = 
(1 — {az))/2 for the z-component of a large spin in the semi- 
classical model, with various initial conditions and e = 1.35. 



We emphasize that a large spin will display the classi- 
cal dynamics of a magnetic moment precessing in a time- 
dependent magnetic field. Plots of the probability as a 
function of time for the semiclassical equations of motion 
of a superradiant ensemble of spins are shown in Fig. [8] 
We see multi-stage transitions similar to the quantum 
case. A semiclassical explanation is as follows. Tran- 
sitions occur when the energy separation between spin 
states equals a multiple of the oscillator frequency. These 
occur at the same times = rjv' given by Eq. (24 1 for 



the quantum case when avoided crossings between adia- 
batic energies occur. 

When the oscillator is initially at rest at its equilibrium 
position </j 0, the initial transition occurs at ^ shown 
in Fig. 8a P(t) oscillates about the regular Landau- 
Zener probability Phz — e^'. Subsequent transitions 
occur at t\ and ^2, with the long-time probability much 
different from Plz- The spin dynamics depend strongly 
on the initial state of the oscillator. Fig. 8b shows the 



transition probability for different initial conditions of the 
oscillator with the same amplitude of oscillation as the 
coherent state studied in the fully quantum-mechanical 
model. 




-1.0 




FIG. 9. Time dependence of the rotation angle in the semi- 
classical model, with various initial conditions and e = 1.35. 

The oscillator dynamics show a similar delay in the 
onset of strong oscillations as in the quantum model, 
shown in Fig. [9j We notice a large increase in the am- 
plitude of oscillations at t\ for the cantilever initially at 
rest at equilibrium. For the cantilever initially oscillat- 
ing with amplitude (/Sq, the behavior is very similar to the 
quantum model with the oscillator initially in a coherent 
state. We observe normal oscillations up to t_i at which 
the amplitude decreases, then a large increase at ii, with 
subequent changes at t2^ t-i- While there is not a strong 
dependence on the initial conditions with the same initial 
energy, there is some variation in maximum amplitude. 



VI. DISCUSSION AND CONCLUSIONS 

We have studied the Landau-Zener dynamics of a tun- 
neling spin rigidly coupled to a torsional oscillator. Start- 
ing with a quantum model describing the low energy dy- 
namics of a tunneling macrospin, we numerically solve 
the time-dependent Schrodinger equation to obtain the 
dynamics of the expectation values of the spin and oscil- 
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lator. We find that when the oscillator is initially in its 
quantum ground state, there are a series of plateaus in 
the staying probability as a function of time. We analyt- 
ically obtain exact probabilities in terms of tunnel split- 
tings of the spin which are dressed by the quantum states 
of the torsional oscillator. These results perfectly fit the 
plateaus obtained from numerical simulations. The oscil- 
lator dynamics show abrupt changes in amplitude which 
occur at the same times as the steps between steps of the 
staying probability. For an oscillator initially in a coher- 
ent state we also find a stepwise staying probability curve, 
but these deviate from the analytical results found for the 
initial ground state because there are multiple occupied 
states of the resonator. The oscillator dynamics continue 
to show changes in amplitude which coincide with the 
steps. We also consider a large number of spins, N, on 
a single oscillator, and find a superradiant enhancement 
of the spin-oscillator coupling which scales as V N. As in 
the Dicke model, the ensemble of spins acts as a single 
large spin. This justifies decoupling quantum averages of 
separate observables in the Heisenberg equations of mo- 
tion, giving semiclassical equations of motion for a large 
spin in a time-dependent effective field which depends on 
the motion of the cantilever. The cantilever experiences 
a harmonic restoring torque but also a driving torque due 
to the dynamics of the large spin. We numerically solve 
the set of coupled equations and compare the results to 
the Schrodinger picture. The spin dynamics show sensi- 
tivity to the initial state of the resonator, although the 
oscillator dynamics are fairly insensitive to this. 

It is important to distinguish the interpretation of 
these results in the context of the system being mea- 
sured. Consider the single molecule magnet grafted to 
a carbon nanotube, depicted in Fig. [Ta| With the sys- 
tem prepared in the spin down state by a strong mag- 
netic field along the negative 2;-axis, the magnetic field 
is swept. If the oscillator is initially in the zero phonon 
state, the first crossing of an occupied energy level with 
an unoccupied level occurs at to = between Eqi and 
Eq-^. Pq = Poo = e~'^°° is the probability that the spin 
will remain in the down state. If the spin remains in the 
down state after the first crossing, it will encounter a sec- 
ond crossing between i^o^ and Ei^ at ti, at which it will 
remain spin-down with probability Poi = e~'^°^ . The to- 
tal probability of the spin remaining spin-down after ti is 
Pi = Poo-foi = e"^'"""'"'^''^''. If the spin reverses at any tk 
it will see no more crossings. When the spin does reverse 
it will exert a torque on the carbon nanotube, exciting a 
phonon mode. The onset of oscillations shows that the 
spin has tunnelled. This provides a method of detecting 
the mechanical quantum state of the nanotube. 

This situation is similar to the recent demonstration 
of electronic readout of nuclear spin states of a terbium- 
based single magnetic molecul^^Zl, Terbium nuclear spin 
3/2 has four possible projections onto the quantization 
axis, each projection providing a different hyperfinc shift 
of the resonance of the Landau-Zener transition of the 
spin of the molecule. Time-resolved measurements show 



an increase in the differential conductance at the time 
the spin makes a transition. This occurs at a different 
value of the external field for each sweep, that depends 
on the quantum state of the nuclear spin. In our model 
the role of nuclear spin states is played by the resonator 



states given by Eq. (23). We, therefore, propose a similar 
experiment in which the field sweep is used to read out 
the quantum state of the mechanical resonator. 

When there is a large number of magnetic molecules 



on a cantilever, as in Fig. lb they will act as a single 
large classical spin R(t). This spin not only responds to 
the external field but also to the motion of the cantilever. 
The latter has been treated in Sec. V as a classical oscil- 
lator described by the angle 4>{t). Such treatment is the 
classical limit of the quantum-mechanical consideration 
in which the cantilever is described by the coherent state, 
Eq. ( 33 1 . Mechanical rotation at an angular frequency (f> 
is equivalent to a magnetic field B^s = 0/7, where 7 
is the gyromagnetic ratio. In turn, the spin dynamics 
act as a driving torque on the cantilever, resulting in 
coupled dynamics which change at the same moments of 
time, tfc, as in the quantum case. The non-linear coupled 
equations of motion lead to the excitations of harmonics 
of the cantilever that correspond to its quantum modes 
in the dynamics described by the Schrodinger equation. 
Higher harmonics are excited with smaller amplitude. 

To put some of these statements into perspective, con- 
sider a spin-10 single molecule magnet grafted to a car- 
bon nanotubdl^n! xhe moment of inertia of the magnetic 
molecule is of the order ~ 10""*^ kg-m^. With a carbon 
nanotube torsional stiffness of fc ^ 10^^^ N-m the sim- 
ple harmonic model gives uJr ^ 1000 GHz, which means 
coupling on the order of A ^ 10"^. Typical phonon fre- 
quencies of carbon nanotubes in the 10-100 GHz range 
would increase the coupling by an order of magnitude. 
Recent observation^ of strong spin-phonon coupling in 
such a system estimates A ~ 0.5. While this is certainly 
large enough to observe the influence of the oscillator on 
spin dynamics, there is no way to directly observe oscil- 
lations in a carbon nanotube. 

If the same spin-10 magnetic molecule were mounted 
on a paddle-shaped torsional resonator of size 20 x 20 x 10 
nm'^ supported by a single carbon nanotube with tor- 

The moment of inertia is 

10~"^^ kg-m'^, which gives 

'^r — \/k/Iz 10^ s~^. The coupling parameter A is 
then on the order of 10^^, which would be too small to 
observe an effect on the spin dynamics. With A/h <C 10^ 
s~^ there should be a detectable delay between the t = 
crossing and the onset of maximal oscillation amplitude. 
With A/h > 10^ s-\ the delay will be undetectable. The 
tunnel splitting can be tuned by orders of magnitude by 
applying a transverse magnetic field. 

A macroscopic resonator in which even small ampli- 
tude oscillations could be observed comes at the expense 
of weak coupling with no observable effect on the spin 
dynamics. In terms of the moment of inertia and tor- 
sional stiffness, the coupling goes as A cx 1/ kl^ , so a 



sional rigidity fc = 10 N-m. 
dominated by the paddle, Iz ^ 
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very small torsional stiffness of fc ~ 10~^^ N-m would 
be needed. One way to overcome this limitation is to 
put a large number of spins on a torsional resonator 
or microcantilever. For a cantilever with dimensions 
1000 X 200 X 100 nm^ we would expect w - 1 GHz with 
Q ^ 500. Single molecule magnets have a diameter on 
the order of 1 nm. It would be possible to place hundreds 
of single molecule magnets on the tip of a nanocantilever 
separated by over 10 nm from their nearest neighbors to 
weaken dipolar interactions. They would act as a single 
large spin due to the collective quantum effect of super- 
radiance. This would increase the coupling by at least an 



order of magnitude, as \sr oc \/N . Therefore it would be 
possible to directly observe the coupled dynamics of the 
magnetization and oscillatory motion in a Landau-Zener 
experiment. 
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